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iNTRorucnoN 

This  report  describes  a  computational  program,  now  in  existence,  for 
estimating  the  accuracy  vith  -idiich  satellite  positions  can  be  calculated 
at  any  time  on  the  basis  of  observations  at  various  stations.  The 
error  in  computed  position  is  due  to  two  main  causes,  at  least  if  sta¬ 
tion  location  uncertainties  are  ignored.  These  are:  (a)  observation 
errors,  (b)  fluctuations  in  the  orbital  decay  rate  due  to  drag.  It  is 
assinaed  here  that  xincertainties  in  the  earth’s  gravitational  field  have 
negligible  effect  on  the  satellite  position. 

The  treatment  of  (A)  follows  standard  statistical  practice:  Orbital 
parameters  are  adjusted  to  a  least— squares  fit  to  the  obseivatlons. 

The  covariances  of  the  estimates  of  the  orbit  parameters  are  then  ob¬ 
tainable  quite  easily  from  the  variances  of  the  observation  errors, 
assuming  that  the  different  meeisurements  are  statistically  independent 
with  zero  bias. 

The  orbital,  parameters  are  seven  in  number,  rather  than  six:  The  (mean) 
orbital  decay  rate  is  a  seventh  parameter  (but  numbered  5  in  §  2)  to  be 
deduced  from  the  obseivations  rather  then  from  prior  knowledge  of  the 
atmosphere.  In  the  case  of  a  received  Doppler  signal,  the  emitting 
frequency,  assxaned  constant  during  a  single  "pass,"  (i.e.,  during  the 
reception  of  the  signal  at  a  single  station  on  a  single  revolution) 
constitutes  an  additional  parameter  per  Doppler  pass.  These  additional 
parameters,  however,  are  easily  eliminated  from  the  calculation,  as  in 
reference  [2]. 


-1- 


The  treatment  of  (b)  is  as  follows:  Fluctuations  in  ’h^--  ie 1  >n 
due  to  drag  are  considered  to  constitute  a  staticiiary  time  s»ri- 1>  wit’n 
exponentially  decreasing  auto-correlation.  The  effect  of  these  irao 
fluctuations  on  altiaide  is  ignored:  Only  the  angulai*  position  aj-ong 
tne  orbit  is  considered  to  be  affected.  A  typical  correlation  time 
lor  the  d -ag  fluctuations  might  be  3  hours.  This  would  correspond  in 
§  6  to  ^  hn  .  The  r.m.s.  value  of  the  drag  fluctuations  may  be 
asstuned  to  be  a  small  fraction,  e.g.,  one  fifteenth,  of  the  drag  at 
perigee.  A  3“0  fluctuation  of  one  fifth  in  density  is  not  unreasonable 
over  short  intervals  (i.e.,  in  a  single  day).  This  is  not  to  be  con¬ 
fused  with  significant  long-term  changes  in  drag  due  to  a  gradual  change 
in  the  relative  positions  of  perigee  and  the  sun. 


§  1.  I£Aar  SQUARES  PROCEDURE 

The  ’'maximum  likelihood"  estimation  of  orbital  parameters  from  a  set  of 
observations  vith  independent  errors  with  various  standard  devia¬ 
tions  requires  the  minimization  of  the  sum  of  the  weighted  squared 

residuals: 

)-q 

(X.1)  tT" 


when  q(?^,t^)  is  the  theoretical  value  for  the  f-th  obseivation,  at 
time  t.  ,  based  on  a  set  X  of  orbital  parameters  X  .  The  index 
<  is  understood  to  run,  say  chronologically,  over  all  observations, 
including  simultaneous  measurements  of  different  physical  quantities 
such  as  range,  azimuth  and  elevation,  in  which  case  thr^e  (successive) 
t^'s  would  be  identical.  The  range  measurements  from 

similar  radar  equipment  will  be  equal,  say  a„  ;  the  o  's  fer  elevation 
measurements  will  all  be  equal,  say  ;  etc. 
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By  linearizing  the  <l(A,t^)'s  as  functions  of  in  x,he  neighborhood 


of  cert^ain  approximate  values  X^  : 


(1.2) 


<l(?v,t^) 


^  ^(X^t  ) 

q(X  ,t^)  +  E  ~  ' 

a 


the  estimation  reduces  to  solution  of  the  following  set  of  simultaneous 
linear  equations: 


(1.3) 


vhere 


2  (K-K)  =  V  , 


& 


r 


(l.U) 


=  Z 


00  '  0  0"  a 


^  dq(X°,t^)  ^(x°,t^) 


00 


a 


BX 


V  =  z 


V. 


a 


1  ^(^*"7^^^) 


dX 


o 


[5,  - 


a 


What  is  of  Interest  to  us  here  is  the  statistical  theorem  that  the 

variances  and  covariances  of  our  estimates  X  of  the  parameters  X 

o  a 

are  just  the  elements  of  the  inverse  of  the  ” information  matrix" 

«c#- 

It  follows  that  if  x  =  x.(X,t*)  is  a  position  coordinate  x,  at 

M.  ^  A 

some  time  t*  ,  computed  from  our  estimated  orbital  parameters  X^  , 


then: 

(1.5) 


\/a  0  ^ 


J 


3x  (X,t*)  ^^(X, t*) 

— — 


0 
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§  evaluation  of  the 


j 

We  proceed  next  to  the  forms  for  appropriate  to  a  particular 

aq(A°,tj) 


a 


will  la<er  b** 


parametrization  of  the  orbit.  The  quantities 

hx. 

related  to  the  's  at  time  by  some  rather  obvious  tri;;onometry, 


and 


O 

will  be  thus  obtainable. 


We  describe  an  orbit,  including  its  decay  due  to  drag,  by  seven 
parameters: 


{?..!) 


ils 

Xg  =  e  cos  p  ,  ^3  =  e  sin  0  t  , 


da 

de 


X^  =  -  fl  sin 


X^  =  i 


Here  a  is  the  initial  semi-major  axis,  a„  an  a-priori  "nominal" 

-I'l 

value  for  a  ;  €  is  the  eccentricity,  p  the  "argument  of  perigee," 

tfj  the  time  st  the  first  ascending  node;  u  is  the  universal  gravita- 
tional  constant  times  the  cir-ss  cf  the  earth;  -  ■—  is  the  rate  of  decay 

Qu 

of  the  semi-major  axis  per  angular  distance  0  from  the  first  ascending 
node;  n  is  the  right-ascension  of  the  ascending  node,  {  the  orbital 
inclination  (i.e.,  the  angle  between  the  orbital  plane  sind  the  equatorial 
plane  taken  as  acute  if  and  only  if  the  orbit  is  Eastward);  is  a 

nominal  orbital  Inclination. 


The  earth's  oblateness,  of  course,  causes  slow  changes  in  the  parameters 
fl  and  p  .  The  amounts  of  their  changes,  however,  are  known  functions 
of  the  other  orbital  parameters  as  well  as  of  the  earth's  oblateness 
coefficient  J  .  which  is  known  with  reasonable  accuracy.  The  accuracy 
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problem  is  therefore  not  appreciably  worsened  by  the  earth's  oblateness 
which  we  shall  therefore  ignore. 

At  each  point  of  a  nominal  orbit  a  right-handed  coordinate  system  is 
used:  x  =  horizontal  In  the  nominal  orbit  plane,  forward; 
y  s=  (horizontal)  perpendicular  to  the  nominal  orbit  plane,  to  the  left; 
z  =  vertically  upward.  It  is  not  difficult  to  see  that  if  denotes 

the  difference  —  (X^)jj  between  actual  and  nominal  parameter  values, 
and  Ay  the  deviation  p>erpendicular  to  the  nominal  plane,  then  (neg¬ 
lecting  second— order  differences); 


^  =  cos  0  AX-  +  sin  6  AX., 

J*  O  • 


r  being  the  local  (radial)  distance  frcaa  the  earth's  center.  Assmaing 
that  the  eccentricity  €  is  small,  we  may  replace  this  by; 


(2.2) 


^  =  cos  6  AX-  +  sin  0  AX- 


Next,  the  radial  dlstemce  r  is  known  as  a  function  of  9  : 


r  = - — ^7^ - r  3f  a(l  -  X  cos  6  -  X  sin  0) 

l  +  €  cos  (0-p)  ''  2  3 


so  that 


=  aj^(l  -  Xj  -  X^  cos  e  -  Xg  sin  e)  , 


(2.5) 


—  =  AX,  -  cos  0  aX_  -  sin  0  AX  , 


ignoring  for  the  moment  the  effect  of  drag. 
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Replacing  a  ry  -r  ‘'j^)  icncrinc  t‘ 


3/2 


(2 


/ 


so  that,  the  Torvarh  ieviaxirn  alcr.r  the 

ry: 

1/? 


(3.1.)  ^=- 


=  3/2 

“>T 


£it  =  -  4  G  -v^.  + 

c: 


+  2  (i  - 


:s  ='  V 


I 


2 


LMSD“70l*013 


§  3.  GALCUIATION  OF  THE  ^  ‘s 

We  come  next  to  the  calculation  of  the  ’s.  In  the  case  of  radar 

range  and  single  measurements  the  qt’s  are  functions  of  the  path  coordi¬ 
nates  X  ,  y  ,  z  ,  as  veil  as  of  the  position  of  the  receiving  station 
relative  to  the  satellite’s  position.  In  the  case  of  Doppler  frequency 
measurement,  the  q  is  a  -^^.ction  of  both  position  and  velocity  along 
the  path  as  well  as  of  the  re'^otive  position  of  the  receiving  station. 
We  shall  postpone  the  discussion  of  Doppler  measurements  until  §  k. 


It  is  clear  that,  for  low  altitude  satellites,  the  partial  derivatives 
vary  much  more  rapidly  with  time  during  a  pass  than  do  the  coordi¬ 


nate  partial  derivatives  x^  ,  y^  ,  z^  themselves.  We  may  thus  get  an 
approximate  set  of  matrix  elements  arising  from  radar  data  by 

writing: 


(5.1) 


^  ^q  ^  .  3q  ''y  .  ^q  hz 

a  a  a  a 


a 


y. 


a 


z 

a 


) 


and  then  replacing  throughout  a  pass  by  its  \ralue  at  the  middle  of 

the  pass.  This  is  equivalent  to  smoothing  alD.  the  radar  data  to  yield 
only  a  single  point  per  pass,  and  relying  on  three  or  more  passes  at 
sensibly  different  points  along  the  orbit  for  the  determination  of  all 
seven  orbital  parameters. 


In  the  case  of  measurements  of  range  S  we  have:  Vq  =  VS  =  u^  ,  a  unit 
vector  along  the  line  of  sight.  Ihe  contribution  to  of  range 

measurements  during  a  single  pass,  assuming  that  the  measurements  have 
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Independent  er^-ors  with  r.m.s.  value  a_  is  thus: 

O 

(3.2)  H^’  =  =  (^)  ^a- S 


s 


when  Q  is  the  matrix  (dyadic)  S  u,  u.  .  Now  if  ♦  denotes  the 

O  ^ 

angular  distance  along  the  orbit,  measured  forward,  from  the  point 
nearest  the  station,  and  if  p  denotes  the  minimum  slant-range  divided 
by  the  earth* s  radius  R  ,  and  if  a  denotes  the  angular  distance  (at 
the  center  of  the  earth)  between  the  station  and  the  orbit  plane,  the 
components  of  the  unit  vector  u^^  are  approximately: 


(3.3) 


-  ,  m  = 


p^  +  cp^ 


p^  +  cp^ 


the  +  sign  to  be  taken  whenever  the  station  lies  to  the  right  of  the 
orbit  (relative  to  the  satellite’s  motion). 

p^rn 

Next,  the  sum  Z  (  )  is  replaced  by  n  /  (  )  dtp  ,  where  n  is  the 

number  of  Independent  radar  measurements  per  geocentric  radian  of  orbit, 
and  is  one  half  of  the  geocentric  angular  interval  under  radar 

observation  from  this  station.  Hence: 


(3.^)  Qc 


„  f  sfJSL 

J 


r  ^  a-  dc 

V 


o%n  rj- 
fi  r  c^'yp  ■ 

J  p^  + 


p^7n  f-^  f^rn  ^ 

r  .J  ^ 

p^  +  cp^  'J  ^  P^  +  ' 
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where  the  off-diagonal  terms  and  vanish  since  th'^* 

inteinral  of  integration  is  symmetric  about  =  0  and  these  integrands 
are  odd  functions  of  4>  . 


Evaluating  the  integrals: 


2n(4.^^-pij0  0 

0 

1 

(5.5) 

fO 

CO 

II 

0  2«  tl' 

P 

±  2n  ^  t 

0  ±  2n  ~  4' 

2n  1  P-^  h'' 

P  V 

V  P  / 

1 

where 


(5.6) 


t  =  Tan 


-  4> 

-1  m 


Turning  next  to  the  angle  measurement^  ve  shall  assume  that  the  r.r„s. 
angular  error  in  the  vertical  plane  is  equal  to  that  in  the  plane 
through  the  line  of  sight  perpendicular  to  the  vertical  plane,  i»e., 
(cos  E)  o,  =  ,  where  E  now  denotes  ele%'ation  angle  and  A  asi- 

muth.  This  assiunption  is  not  valid  for  very  low  elevations  where 
is  substantially  greater  th.an  .  However,  these  low  elevation 
measurements  are  usuailv  excluded:  i.e.,  the  value  of  ^  corresi>onds 
to  the  miaimum  allcvabie  E  which  is  greater  than  zero. 


Introducing  two  more  unit  vectors  and 

a  perpendicular  tri ' 1,  the  contribution  of 


u  ,  which  with  u,  form 

3  1 

angular  errors  to  M  _  is 


A 


-ic- 
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when  S  ,  the  instantaneous  slant-i-ange,  Is  approximately  R  -yj  +(p^ 
thus  the  angle  contribution  to  M  „  may  be  written: 


(5-7) 


M 


(E) 


=  ( M 

VeJ 


C  *  S  * 


where 


(Ug  Ug  -»•  U3  U3)&P 


P^  +  <P‘ 


>« 


Introducing  the  components  and  >  ^3  t  ^3  ^2 

and  Ug  and  making  use  of  the  identities  ^  1  ,  etc.,  and 

^1^1  +  4 ”^2  +  -^3 f  etc.,  we  obtain 


.8)  ej,= 


i_,2 


-<p. 


m 


J  +cp^ 

~<P 

m 


r”"  .  r  ""  -”1  -i 

^  '”7  p2+<p= 

-q)  -<!) 

TK 


dcp 


cp 

--J  -"7 

“'P-m 

-/ 


m 

n  /  — z — ^  dcp 


-<P 


p2+q)2 


m  n. 


dcp 


-^1 


-CP 

m 


TP, 


n 


?n 


P^+cp® 


(icp 


P^+cp2 


dcp 


m  1-n^ 
_ 1 

p2+cp2 


dq) 


sin  2  \pj 


0 


0 


0 


2ni 

P 


(..j 


sin  2  i'  T 


0 


7U7 


-iL -  „  _  s!„  2j 

S2i  -  7  +  i  sl„  2  i 
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The  total  matrix  M  _  from  all  radar  observations  is  now: 

03 


(3.9) 


where  the  index  k  runs  over  the  set  of  passes  and  where  the  parameters 
<p^  ,  a  and  p  associated  with  each  pass  (each  k  )  may  be  deduced  from 
the  three  quantities:  Local  altitude  h^  ,  maximum  elevation  ,  and 
minimum  allowable  elevation  E  . 


Thus 


(3.10) 


(3.11) 


where 


and  so 
(3.12) 


<p  c  cos"^ 


cos  g* 
cos  a 


y 


COS  a’  = 


R  cos^  E  +  sin  E  -v/(R+hJ.-  - 

o  o  V  ^  k' _ 


cos^  E 


R  +  h. 


R 


CD  =  cos 
m 


-1 


cos^  E  +  sin  E  "v/(R+h.)^  -  cos^ 
o  o  V  « 


E, 


R  cos^  E^  +  sin  E, 


R^  cos^  E, 
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§  U.  THE  INCUJSION  OF  DOPPLER  DATA 

A  single  Doppler  pass  for  a  low-altitude  satellite  gives,  essentiall;/, 
a  measurement  only  of  the  three  quantities:  Minimtom  slant— range  , 
time  of  minimum  slant— range,  and  speed  . 

It  is  shown  in  reference  [2]  that  if  frequency  f  is  measured  with 
standard  deviation  and  if  there  is  negligible  drift  in  transmitted 

frequency  during  a  single  pass,  but  if  the  actual  transmitted  f  is 
treated  as  an  imknown  constant,  the  information  matrix  relative  to  the 


To  convert  this  into  an  information  matrijc  relative  to  our  seven 

parameters  we  denote  the  maximum  elevation  during  pass  k  by  E,  . 

ct  '  ■  *  fe  ■ 

and  observe  that: 


(4.2) 


iS,  =  dz  sin  E.  ± 
k  K 


dv  cos  E. 


according  as  the  station  lies  to  the 


fS'n’l 


the  orbii 


-1-^- 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  orsd  SPACE  Div'iSiCN 


I 


TMCT^'7/%1  ./ni  ^ 
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Next  we  Introduce  the  notation: 


(^.5) 

so  that,  since 


,a 

Uy 


V  S  re 


~  (a  )  ^  (0-  3)J 


■(A) 


1/2 


^  ~  2  ■*■  ^2  0  +  sin  e 


ve  have: 


r 


(•(.t) 


■< 


Vj,  at  COS  e 


v^  a  Sin  e 


%  =  0 


V5  =9 


Vg  =  v^  a  0 


V 
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Making  use  of  (4.2),  (4.9),  and  the  notation  of  (2.5),  and  replacing 
ajj  by  R  ,  the  sHtrlx  convertible  into  the  7x7  matrix: 


M 


(±  \  sin  .  (±  cos  +  z 


"  a  k 
1/2 


3 


+  2R 


(S)  «[■.(. 


^  y  cos  E.  +  z  sin  E^ 
P  «  p  k 


+  V, 


cos  E. +z  sin  E^ 
k  a  k 


) 

a 


.  ji  JD)  ^  ^  ^  r3  Jd) 

R  33  a  p  p  a  p 


Replacing  v  by  ,  this  is: 


%  ' 


c^R 


='3  *  Bln  ^ 


yg  Bln  Ej 


H  "a  ^  -  h 


±  yp  BOB  +  Bg  nln  Ej) 


■*■  y-v  E,  +  z  sin  E- 

p  \  a  e  a  ife. 


where 


(4.6)  -< 


K  ~  +  ■!  si” 

=  ii^(^  -|=ln 

H  X  '  “  ~ 

^  ~1  2^  -  ^  sin  4 
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§5.  THE  EFFECT  OF  RANDOM  FLUCTUATIONS  IN  THE  DRAG 


Suppose  that  random  fluctuations  In  the  air  drag  have  the  effect  of 
perturbing  the  time  of  arrival  at  a  given  orbital  position  vithout 
significantly  perturbing  the  altitude. 


A  measurement  is  now  related  to  a  set  of  parameters 

(5.1)  =  q{^  ,t^)  + — — 5x{t^)  +  n. 


\ 

a 


by: 


9 


where  is  a  standard  normal  variable,  and  where  Bx(t^)  is  the  per¬ 

turbed  position  along  the  orbit  at  time  t^  .  The  6x(t^)*8  corresponding 
to  different  times  t^  are  correlated  through  their  dependence  on 
previous  air  drag  fluctuations  (see  below). 


The  least-squares  fitting  procedure  described  in  §  1,  which  tacitly 
Ignores  any  drag  fluctuations,  leads  (as  we  have  seen)  to 


(5.2) 


A 


(M- 


Linearization  of  (5.T)  neighborhood  of  A  =  X®  and  substitution 

into  (5.2)  leads  to; 


(5.5) 


P  Op  i 


Sq(X°,t^ 


aq(X°,t^) 


*  5x 


6x(t^)  + 
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But  L  — 


1  dq()>°,tj) 


Equation  (5-5)  thus  simplifies 


(5.i»)  X  -X  =Z(M'^)_Z—  - 

a  a  Q  ^  i  °4 

p  3 


K-  +  — 

i  o. 


1 


6x(tjl 


Denoting  the  covariance  of  the  quantities  ox(t  )  ,  6x(t.)  by  u 

1  ^  i  j 

and  assuming  that  these  quantities  are  noraally  distributed  vith  mean 

zero  and,  of  course,  independently  of  the  measurement  errors  -• 

foUovs  that  the  covariance  of  the  estimates  X  is; 

a 


1  3q(^°,t^) 


Sq(X°,tp  hjjt.) 

£  i  Oj  £-£  5x  Sx 


This  simnlifies  to 


(5-5) 


^  )(>'o 

a  p 


-  v]  ■ 


^  oS  ''  ^00 


■where 


(5.6) 


E  E  U. 

f 


^2  ^g(tp^ 


0\  erf 


Kext,  as  in  §  5.  group  the  summation  over  observations  into  sumnations 

'hr 

over  Dasses.  Like  -^r-  ,  u,,  mav  be  reftardei  as  effectivelv  constant 

<y\  tj 

vhen  i  irtns  over  obser’/ations  during  a  sinrJle  pass  and  j  runs  over 
observations  during  another  pass,  at  least  if  the  correlation  time  of  the 
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6x(t)*8  is  fairly  long,  as  ve  may  expect.  Using  the  particular  form 
of  (4.5)  and  the  matrices  (3-^)  and  (3‘‘^)  vherein  x  is  not  coupled 
to  y  ,  z  ,  the  contribution  of  a  single  pass  k  xo  the  sum 

n  ai(x‘^,t  )  d<l(t  ) 

^  a 

is  Just 

1  .  1  . 

— r  A.  - - —  =  —  A.  X 

K  “ 

where 


(5.7) 


Equation  (5*6)  thus  siii5)lifies  to 


(5.8) 


N 


OfP 


Z  Z 


kg 


A.  x„  x^ 


*  £  a; 


2  "1  a  ^ 


N 


where  the  sinoEitlcxis  E  and  Z  are  now  suaanations  over  different  nasses, 

k  t 

Before  evaluating  the  covariances  of  the  6x(tj)*s  ve  must  note 

that  the  true  position  is: 


(5.9) 


x(t)  »  x(X,t)  +  6x(t) 


where  X  are  true  parameters  at  scms  time  t  .  and  6x(t)  is  the 

O' 

accxifflulated  position  error  since  t  ,  expressible  in  the  form  (cf.  §6): 


(5.10) 


6x(t)  =r  /  (t-t*)  6x(t»)  dt»  . 

"'t 
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The  error  in  position  (Jetennination  along  the  orbit  at  tinje  t*  is: 
(5.11)  €^(t*)  =  x(A,t*)  -  x(A,t*)  -  6x(t*) 


Choosing  t^  to  coincide  with  t*  ,  we  have  6x(t*)  =  0  and 


(5.12) 


x(t*) 


a  3 


(M-*) 

ap  03 


^(A,t*)  ^(A,t*) 


^A 


a 


3 


The  first  term  in  the  [  3  i-eproduces  the  variance  due  to  measurement 
errors  described  in  §  1.  The  second  term  yields  the  contribution  of 
drag  fluctuations. 


§  6.  evaluation  of  the 


It  is  easy  to  show  from  energy  considerations  that  a  small  change  in 
velocity  along  a  near— circular  orbit  yields  a  change  in  semi-major 

axis  given  by: 


(6.1) 


If  the  change 


1  dr 
r  and  —  .yr- 

r  de 

notation  of  S 


occurs  at  angular  position  A’  from  the  equator,  since 
must  renain  unchanged  at  this  position  ve  have,  in  the 
2: 


(6.2) 


CCS  e’  -f  iiA  sin  e'  = 


sin  8*  +  cos  ?’  =  C  , 


—1 
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80  that 


(6.5) 


M  a  2  —  cos  e » 
2  V 

AXg  =  2  ^  sin  e» 


Substitution  into  (2.4)  yields  a  change  in  subsequent  distance  along 
the  orbit,  given  by: 

(6.1,)  =  ^ 


-  3(0  ~  0’)  +  4  sin  (e  -  e’) 


Now  the  reteirdatlon  due  to  drag  is  given  by 


(6.5) 


dv  „  2 

^=-Bpv2  , 


where  p  is  air~density  and  B  *  Cj^a/w  ,  being  the  satellite 

drag-coefficient,  A  the  reference  area,  and  m  the  satellite  mass. 


The  effect  on  position  x(8)  due  to  drag  between  0*  and  0  is  thus 
expressible  by; 


AxCo) 

a 


9  p 

=  /  Bpv  j  5(9  -  0* )  -  U  sin  (0  -  0'  )  dt*  , 

“  i=9.  L  J 


V  de*  V  ... 

where  -rrr  =  —  »  so  that 

at- 


(6.6) 


Ax(9) 


0 

n 


J  (BP 

V 


;-(e  -  e»)  -  4  sin  {0  -  6») 


de* 
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(6.9) 


The  effect  of  an  exponentially  decreading  density-altitude  relation¬ 
ship  is  accounted  for  by  the  parameter  of  §  2.  The  effect  of  a 

deviation  6p(e*)  in  density,  at  angular  position  0*  ,  from  that 
given  by  the  exponential  formula  yields  a  change  in  position  6x(0) 
from  that  corresponding  to  the  set  A  of  parameters  vhich  fit  per¬ 
fectly  at  0*  ,  namely: 


(6.7) 


6x(e) 


B  bp(0» ) 


5  (9  -  e' )  -  4  bxn  (0  -  0» ) 


d0*  . 


Next,  vs  shall  assume  that  6p(0’)  is  a  stationary  time-series  with 
negative  exponential  auto-correlation,  so  that 

(6.8)  S'|^5p(0i)  MQj,) 

vhere  is  the  variance  of  8p  and  0  the  non-dimensional 

P  c 

’’correlation  time’*  of  the  time  series. 


} 


s=  e 
P 


- 101-92 1 /Oe 


Assuming  that  0^  >  2n  ,  say,  and  that  the  sum  (5*8)  contains  passes 
over  several  revolutions,  the  second  terra  in  the  square  bracket  Inside 
the  integral  (6.7)  is  unimportamt  by  comparison  with  the  first  term. 
Omitting  the  unimportant  term  and  using  (6.8)  we  obtain: 


6x(0^)  5x(0^) 


Finally,  the  double  integral  in  (6.9)  may  be  evaluated  to  yield: 

(6.10)  =f 


-  f(|0-»^-oJ)  ±  5(0*- 0^) 
±  5(O*~0j^)  r(!  0*-  9^  I  )  ■ 
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where  the  in  siffn  oovrespond  to  o*  p  =u5d  ^ 

k  t 

respect lve3.y  and  where; 


\ 


H 

I 


— fl/ 0^ 

e'^  e  '  de  ^ 


i  6-" 


:5  0  +  fv  p  of 

c  c 


y 

/ 


-e/8,. 


e""  e  do  =  e. 


e'^  - 


6  9®  +  6  8^  e 
c  c 


2  e  6?  +  2  0? 


~6/e, 


-e/e. 


2  9*-  e 
c 
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